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Abstract . The symmetrical problem of two transverse cracks in an 
elastic strip with reinforced surfaces is formulated in terms of 
a singular integral equation. The special cases of one central 
crack or two edge cracks are discussed. Numerical methods for 
solving the problems with internal cracks are outlined and stress 
intensity factors are presented for various geometries and degrees 
of surface reinforcement. 


1. INTRODUCTION 

In a recent paper [1] Gupta and Erdogan have treated the 
problem of symmetrical transverse' cracks in an elastic strip by 
means of singular integral equations. References to alternative 
approaches to similar problems are given in [1] . In this paper 
we consider transverse cracks in an elastic strip with reinforced 
surfaces. The bending stiffness of the surface reinforcement is 
neglected, and the effect of the reinforcement is accounted for by 
a modification of the boundary conditions for the elastic material 
The problem is formulated as a singular integral equation with the 
derivative of the normal crack surface displacement as the unknown 
function. A numerical solution is carried out for several geom- 
etries and degrees of reinforcement, and the stress intensity 
factors are given. 

* This work was partially supported by a grant from NASA-Langley 
NGR- 39-0 07-0 11 . 

**Present address: Structural Research Laboratory, Technical 

University of Denmark, 2800 Lyngby, Denmark. 
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2. THE INTEGRAL EQUATIONS 

Consider the strip shown in Figure 1. The geometry and the 
external loads are symmetric with respect to x=0 and y=0. The 
sides of the strip are reinforced with membranes with finite 
stretching stiffness T. The bending stiffness of the membranes 
is neglected. By using symmetry considerations it is seen that 
the corresponding perturbation problem can be formulated with the 
following boundary conditions: 


u(0,y) = 0 , 

a xy (0 ' y) = 0 

/ 

0 < 

y < °°; 

(1. a,b) 

o xx (h,y) =• o , 

°xy (h,y) = T 
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— P v (h,y) , 

9 y 
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y < 
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a yy (x,°°) = o , 
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o xy (x,0) = 0 , 
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v (x, 0) = 0 , 
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u and v are the x- and y-components of the displacement vector. 
Following [1] and [2] the solution is given in the form 
u(x,y) = - | 


oo 

[ f ( s ) - g(s) ] sinh(sx) 

O ' 

+ xg (s) cosh (sx) } cos (sy) ds 



m(p) 

P 


- py)e ^ y sin(px)dp 
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v (x,y) 
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+ 


+ py) e py cos (px) dp 

P ^ 

CO 

{i [f(s) + g (s) ] cosh (sx) 

o s ^ 

+ xg (s) sinh (sx) } sin(sy)ds 


(7.a,b) 


When the solution is expressed in this form the field equations 
and the conditions (1) , (3) and (4) are satisfied. The material 
parameter k is' defined by k = 3-4v for plane strain and 
k= (3-v)/(l+v) for plane stress, where v is Poisson's ratio. The 
functions m(p), f(s) and g(s) will be determined from the remaining 
boundary conditions (2) and (5). From (7) the stresses are found 
to be 


2y a xx 


2y a yy 


1 

2y 



+ 


+ 





m(p) (1 - py) e”^cos (px) dp 

[f (s) cosh (sx) + sxg (s) sinh (sx) ] cos (sy) ds , 
m (p) ( 1 + py ) e -p) ^cos (px) dp 

{ [f (s) + 2g (s) ] cosh (sx) 

+ sxg (s) sinh (sx) } cos(sy)ds , 
ypm(p)e p) ^sin(px)dp 

{ [ f ( s ) + g (s) ] sinh (sx) 

+ sxg(s) cosh(sx) } sin(sy)ds (8.a-c) 
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where y is the shear modulus. By transforming (2) and (8) the 
following equations are obtained 


4 s 

f(s)cosh(sh) + shg (s) sinh (sh) = - — 


— 2 COS dp ' 
o (p Z + S Z ) 


f(s)sinh(sh) + g (s) [sinh (sh) + (1 + $) shcosh (sh) ] 


_ 4s f p 2 m(p) 


•n j ,7 
J o (p z 


1 ^ p ^ ■ 2 sin(ph)dp + $sh % - S m ^ p l- cos(ph)dp 
+ s 2 r U o p z + s z 


(9 .a,b) 


where the degree of reinforcement $ is defined by 


$ = T = 

4yh 


T 

2 for plane stress 


(1 — v ^ ) T 

2 - — — - — for plane strain 


(10) 


Eh 


Straightforward use of the mixed boundary condition (5) for 
the determination of the last unknown function m(p) would lead to 
a set of dual integral equations. Here the method from [1] is 
used. A new unknown function- G(x) is defined as 


G(x) = v(x,0) , 


0 < x < h 


( 11 ) 


By use of (5.b) and (7.b) the following relation is obtained 


2 ( b 

m ( p ) = - G (t) sin (pt) dt 


( 12 ) 


Substituting for m(p) by means of (12) in (9) and (8.b) the last 
boundary condition (5. a) takes the form 


fb 

H(x,t)G(t)dt = - Tra (x) , a< x< b (13) 

•'a ^ 

Noting that the first term of (8.b) gives rise to a Cauchy- type 
singularity, (13) may be expressed as 


4 ‘ 
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J [ t^x + tfe + k(x,t) 1 G(t)dt = -- 7Ta(x) , 

a < x < b , (14) 

where 

f°° 

k(x,t) = [K(x,t,s)e tS - K(x,-t,s)e ts ]ds 
■*0 

K(x,t,s) = { [ (1 - $ sh) (sh - 2) - (2sh + $sh - 3) s (h-t) 

- ( (1 - $sh) (sh + 2) - (1 - $sh) s (h-t) )e 2sb ] cosh (sx) 

- [1 - $sh - 2s (h-t) + (1 - $sh) e -2sh ]sxsinh (sx) } 

/ [sinh(2sh) + 2sh + 2$shcosh 2 (sh) ] (15. a, b) 

The index of the singular integral equation is +1, and the neces- 
sary extra condition is given by (6) , 

f b 

G(x)dx =0 (16) 

a 

For 0<a<b<h the kernel k(X/t) is bounded. Thus (14) is an 
ordinary singular integral equation with the fundamental function 

w (x) = [(X- a) (b — x) ] — 1/2 (17) 

A numerical solution can then be obtained by the method described 
in [3] and [4 ] . 

For a=0, b<h the problem reduces to a single central crack. 
No additional condition is present in this case. The term l/(t+x) 
now becomes unbounded for x,t-*0. This influences the fundamental 
function. k(x,t) is bounded and the fundamental function may be 
determined in a straightforward manner by variable transformation 
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[5] or reflection in the line x=0 [1]*. The result is 


w(x) = (b 2 - x 2 ) 1/2 


(18) 


3 . EDGE CRACKS 

For a>0, b=h the character of the integral equation changes 
as the function k(x,t) becomes unbounded for x,t-*-h. In the fol- 
lowing calculations there is an essential difference between the 
cases $ = 0 and $ /0. This is caused by the denominator of (15. b) , 
where the term of highest order for s -*• 00 is multiplied by $. As 
the case $ =0 has been treated in detail in [1], only $ / 0 will 
be considered here. In (15. b) terms which for x,t-»-h become of 
order less than 1/s for s -*■“ are extracted and integrated analyti- 


cally . 

, This is done by use of the formula [6] 
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^ . J 

Jtl 

* ds = — — 
dt m J 

f" -s (2h-t) 

© i 

o 

sinh (sx) 

> 

cosh (sx) 
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_ d m 1 
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r ■> 
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dt m 

2 2 1 

_(2h - t) - x z 

2h- t 

v J 

t 


= . I [ 1 . + 1 ] 

2 dt m t - (2h-x) t - (2h + x) 1 


(19) 


It appears that only the terms which originally contained the 
factor $ contribute to the singular part of the kernel. Only 
these terms are extracted in the following formulation. 

| ' rh - t - uh-x) + k f (x ' t) J G(t> = - ” (x) - 

cl 

a < x < h , (20) 

*The reflection technique applied in [1] requires a slight modi- 
fication of the kernel, which would otherwise be defined by a 
divergent integral. 
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where 


k f (x,t) t+x t - (2h + x) 


[K(x,t,s)e ts - K(x,-t,s)e ts - (x , t ,s) e ts ] ds 


K^Xft/S) = 2e S ^{ [- sh + 2 - s (h-t) ] cosh (sx) + sx sinh(sx)} 

(21. a, b) 


The first two terms in the kernel of (20) constitute a 
generalized Cauchy kernel. It is noted that this part of the 
kernel — and thereby the form of the singularities — is independent 
of $, provided $ ^0. Because of the simple form of the generalized 
Cauchy kernel the fundamental function can be determined directly 
by change of variables [5] . The result is similar to (18) . 

w(x) = [(a-h) 2 - (x-h) 2 ]“ 1/2 (22) 

In the case where a strip with edge cracks is loaded with 
stresses not acting on the crack surfaces, the problem formulated 
here corresponds to membranes broken at y =0. This is seen in the 
following way. First the boundary value problem for the uncracked 
strip is solved. This solution yields the normal stress a(x) on 
the crack surfaces.. For reasons of continuity the membrane force 
at y = 0 is 

s = ? nr lim a(x) (23) 

M 0 x+h 

When the perturbation problem is solved with the load -c(x) on the 
crack surfaces, continuity requires the membrane to be loaded with 
the force -S. Thus the resulting membrane force is zero at y = 0 
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and no interaction takes place between the membrane for y > 0 and 
y < 0. 


4. NUMERICAL METHOD 

The objective is to determine G(t) from the singular integral 
equation (14) . First we consider the case 0 < a <b < h. The ker- 
nel then has a normal Cauchy singularity and we normalize the 
variables by defining t and E, 

x = (b+a)/2 + £(b-a)/2 

t = (b+a)/2 + t (b-a) /2 (24. a, b) 

Furthermore we introduce the notation 

G(t) = g(T) (1 -t 2 )" 1/2 (25) 

where we have used the fundamental function (17) . (14) may now be 

expressed as 

1 M(S,T)g(T) (1- T 2 )" 1/2 dT = - ^ Trcr(S) , -1<C<1 (26) 

a(£) is used for convenience for a(x) in connection with (24. a). 

M($,T) = (T-O' 1 + (T+El^ + yk'lx.t) (27) 

(16) takes the form 

| g(T) (1- T 2 )" 1/2 dT = 0 (28) 

Based on the assumption that g(x) can be approximated by an (n-1) 
order polynomial we obtain [3] 

I M( ?i ,T )g(T ) .= - no(E i ) 

3=1 J J 
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n 


l g(T.) = o 
j=i 3 


h = 


x 1 , 2 , ... ,n 1 


T . = COS x ~— 

J 2n 


(— Tr) , j = 1 , 2 , . . . ,n 


(29 . a-d) 


(29) gives n linear equations for the determination of the n 
values g(tj). Special importance is attached to the values g(l) 
and g(-l) giving the stress intensity factors in dimensionless 
form. As proved in [4], g(l) is determined by 

g(i> = K l g(T -i ) (30) 

n j=l sin [ 4 “ IT] 3 

g(-l) is found from a similar formula with g(x^) exchanged by 
9 (T n+l-j ) - 

The case a= 0, b< h is solved by a similar method. We now 
define t and £ by 


x = b£ , 


t = bx 


(31. a, b) 


(25) remains valid, and the integral equation (14) becomes 

1 N(£,x)g(x) (1- T 2 )" 1/2 dx = - irc(£) , 0<£<1 (32) 


N (£, x) = (t -£) -1 + (x+£) 1 + bk (x, t) 


(33) 


By using the fact that g(x) and N(£,x) are odd functions of x we 
may use the formulas (29) to evaluate the integral in (32) . We 
use 2n instead of n in the formulas (29) and thereby obtain: 

n 1+ 

l N(£.,x )g(x.) = - ±+K 


j=l 


■3T 2na(? i ) 
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£ ^ cos ( 2 n u ) 9 ^ 1, 2 , . . . ,n 

t j - cos ( tt), j = 1,2, ... ,n (34. a. - c ) 

g(l) is obtained from (30). 

5. RESULTS 

First we consider the case shown in Figure 1. The value of 
a is kept constant and b and 0 are varied. The Tables 1 and 2 
contain the stress intensity factors for constant a. When compar- 
ing results for different values of b it must be kept in mind that 
the stress intensity factors are normalized with respect to dif- 
ferent values of l. For a= 0.2h, b= 0.6h and $= 0 the following 
results are given in [1] : 

k(a)/a/£ = 1.1102 , k (b )/a/l = 1.0961 

Table 3 contains results for a central crack. For $ = 0 the 
results may be compared with similar data given in [8] and [9]. 
With the exception of b= 0.9h there is agreement to within one 
unit of the last digit. For $ ^ 0 curves are given in [8] . Within 
the error of this kind of comparison the agreement is good. 
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Table 1. The stress intensity factor at a for a= 0.2h. 


a/h=0 . 2 


k (a) /a/£ 

, & = 

(b-a)/2 


b/h 

$ = 0 

$ = 0.2 

$ = 1.0 

$ = 2.0 

4>= 1000 

0.3 

1.009 

1.008 

1.007 

1.006 

1.005 

0.4 

1.031 

1.029 

1.024 

1.020 

1.013 

0.5 

1.064 

1.058 

1.046 

1.039 

1.020 

0.6 

1.109 

1.097 

1.073 

. 1.059 

1.025 

0.7 

1.169 

1.148 

1.105 

1.081 

1.025 

0.8 

1.253 

1.213 

1.140 

1.102 

1.017 

0.9 

1.38 

1.30 

1.18 

1.12 

1.00 

Table 2. 

The stress intensity factor 

at b for c 

= 0 . 2h. 



k(b )/a/l 

/ A = 

(b-a)/2 


b/h 

4> = 0 

$ = 0.2 

$ = 1.0 

$ = 2.0 

$= 1000 

0.3 

1.008 

1.007 

1.006 

1.005 

1.004 

0.4 

1.026 

1.024 

1.018 

1.015 

. 1.007 

0.5 

1.053 

1.047 

1.034 

1.026 

1.008 

0.6 

1.094 

1.080 

1.053 

1.037 

1.001 

0.7 

1.156 

1.127 

1.074 

1.046 

0.984 

0.8 

1.263 

1.197 

1.094 

1.045 

0.946 

0.9 

1.51 

1.31 

1.09 

1.00 

0.85 

Table 3, 

Stress intensity factor at 1 

for a = 0. 

a/h = 0 

k(b)/cr/b 

b/h 

$ = 0 

$ = 0.2 

$ = 1.0 

4> = 2.0 

$ = 1000 

0.1 

1.006 

1.005 

1.003 

1.002 

0.988 

0.2 

1.025 

1.021 

1.011 

1.006 

0.992' 

0.3 

1.058 

1.048 

1.026 

1.013 

0.981 

0.4 

1.110 

1.089 

1.047 

1.023 

0.965 

0.5 

1.187 

1.150 

1.075 

1.035 

0.942 

0.6 

1.304 

1.237 

1.112 

1.048 

0.911 

0.7 

1.488 

1.366 

1.156 

1.058 

0.867 

0.8 

1.817 

1.567 

1.201 

1.055 

0.802 

0.9 

2.57 

1. 89 

1.21 

1.00 

0.68 
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